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We present the calculation of the elastic and inelastic high-energy small-angle electron-positron scattering 
with a per mille accuracy. 



1. Introduction 

A considerable attention has been recently de- 
voted to the Bhabha process 0-1 . The reached 
accuracy is however still inadequate with re- 
spect to the experimental one [Q. According to 
these evaluations the theoretical estimates are 
still incomplete; moreover, they are somewhat 
larger (~ a factor 2) than the projected theoret- 
ical and experimental precision Q and are com- 
parable to the currently published experimental 
precision. 

The process that will be considered in this work 
is that of Bhabha scattering when electrons and 
positrons are emitted at small angles with respect 
to the initial electron and positron directions. We 
have examined the radiative processes inclusively 
accompanying the main e+e"^ e"'"e~ reaction at 
high energies, when both the scattered electron 
and positron are tagged within the counter aper- 
ture. 

We assume that the center-of-mass energies 
are within the range of the LEP collider 2e = 
•^s = 90 — 200 GeV and the scattering angles 
are within the range 6 ~ 10 — 150 mrad. We 
assume that the charged-particle detectors have 
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the following polar angle cuts: 

0-01 < 0, < 0.1 rad , (1) 

where Pi, {p2, 'Z2 ) the momenta of the 
initial and of the scattered electron (positron) in 
the center-of-mass frame. 

In this talk we present the results of our 
calculations of the electron-positron scattering 
cross-section with an accuracy of 0(0.1%). The 
squared matrix elements of the various exclusive 
processes inclusively contributing to the e+e~ 
e^e~ reaction are integrated in order to define an 
experimentally measurable cross-section accord- 
ing to suitable restrictions on the angles and en- 
ergies of the detected particles. The different con- 
tributions to the electron and positron distribu- 
tions, needed for the required accuracy, are pre- 
sented using analytical expressions. 

In order to define the angular range of interest 
and the implications on the required accuracy, let 
us first briefly discuss, in a general way, the angle- 
dependent corrections to the cross-section. 

We consider e+e" scattering at angles as de- 
fined in Eq. (0). Within this region, if one ex- 
presses the cross-section by means of a series ex- 
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pansion in terms of angles, the main contribution 
to the cross-section Aa/dO'^ comes from the di- 
agrams for the scattering amplitudes containing 
one exchanged photon in the t-channel. These di- 
agrams, as well known, show a singularity of the 
type for — > 0, e.g. 

d^ ^ ■ 

Let us now estimate the correction of order 9^ 
to this contribution. If 



(2) 



then, after integration over 9^ in the angular 
range of Eq. (nl), we obtain: 



da 
d02 



(3) 



We see that, for 6',-„i„ — 50mrad and 0„,ax = 
150 mrad (we have taken the case where the 9^ 
corrections are maximal), the relative contribu- 
tion of the 9^ terms is about 5x 10~'^ci. Therefore, 
the terms of relative order 9^ must be kept only 
in the Born cross-section where the coefficient ci 
is not small. 

A detailed derivation of these results has 
been reported elsewhere (see Q and references 
therein) . 

2. Born cross section and one-loop virtual 
and soft corrections 

The Born cross-section for Bhabha scattering 
within the Standard Model is well known ||^. In 
the small-angle limit is 



da^ 
9d9 



£2614 



(4) 



where e = \fsl'2- is the electron or positron initial 
energy and the weak correction term (5„oak, con- 
nected with diagrams with Z^'-boson exchange. 

In the pure QED case one-loop radiative cor- 
rections to Bhabha cross-section were calculated 



a long time ago js) . Taking into account a contri- 
bution from soft-photon emission with energy less 
than a given finite threshold Ae, we have here for 
the cross-section rfcTg^^) , in the one-loop approx- 
imation for the limiting case of small scattering 
angles: 



dc 



dc 



(1 



n(t))- 

3 



(l + '5), (5) 



2(1 - L)ln — + -L 
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The photon vacuum polarization function Tl{t) 
(n(s)) is defined as follows: 



^ ' TT V 3 9 



4 TT 



where 



In- 



(6) 



(7) 



and we took into account the leading part of the 
two-loop contribution in the polarization oper- 
ator. In the Standard Model, St contains con- 
tributions of muons, tau-leptons, PF-bosons and 
hadrons. For numerical calculations we use for 
n(i) the results of Rcf . §. 

Taking into account that the large contribution 
proportional to In A disappears when we add the 
cross-section for the hard emission, one can verify 
that terms of relative order 9^ can be neglected. 
Therefore we will omit in higher orders the anni- 
hilation diagrams and multiple-photon exchange 
diagrams in the scattering channel. The sec- 
ond simplification is justified by the generalized 
eikonal representation for small-angle scattering 
amplitudes. In particular, for the case of elastic 
processes we have |lOi : 



Ais,t) = Ao{s,t) Flit) (1 -n(<))-i e^'^W 
'a Q^- ' 



1 + 



S > > TO^ 



(8) 
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where Ao{s, t) is the Born amphtude, Fi{t) is the 
Dirac form factor and (p{t) — —a ln((5^/A^) is 
the Coulomb phase, A is the photon mass aux- 
ihary parameter. We may consider the eikonal 
representation as correct within the required ac- 
curacyQ. 

Let us now introduce the dimensionless quan- 
tity ^ Ql cr„,p/(47rQ;2), with Ql = e'^Bf, where 
CTexp is the Bhabha-process cross-section inte- 
grated over the typical experimental energy and 
angular ranges^: 

^2 J dxi J dx2 <d{xiX2 - Xc) 



X / d^q^ / d^q^ e=2 



(q^ ,X2) e- [q^ ,xi)+X 



(9) 



where Xi.2, qi2 ^I'c the energy fractions and the 
transverse components of the momenta of the 
electron and positron in the final state, sxc is 
the experimental cut-off on their invariant mass 
squared and the functions do take into account 
the angular cuts (|l]): 



= 6(03- 



e§ = 6(614 - 



) e(^-0i), 
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In the case of a symmetrical angular acceptance 
(we restrict ourselves further to this case only) we 
have: 



?3, 



(11) 



We will present E as the sum of various contri- 
butions: 

= ^ooil + So + S^ + S'''' + 6'-^''' 

+ 6'-' + 6^^^'-'), Soo = l-p-^ (12) 

Result obtained in paper wc believe, is incorrect. It 
contradicts to the well estaoushed result of D. Yennie et 
al. [hd about cancelation of infrared singularities. 
^Really this quantity corresponds to some ideal detectors. 
It is intended for comparisons with the results of Monte 
Carlo event generators. 



where Eq stands for a modified Born contribution, 
E'' for a contribution of one-photon emission (real 
and virtual) and so on. The values of the 5' as 
function of Xc are given in Table 1 (see below). 
We shall slightly modify the perturbation theory, 
using the full propagator for the t-channel pho- 
ton, which takes into account the growth of the 
electric charge at small distances. By integrating 
with this convention, we obtain: 

el 

^a=Olj ^(1 - n(t))-2 + Eh^ + Ee, (13) 
ei 

where Yiw is the correction due to the weak in- 
teraction: 



^2 



d^ 

614 



(5woak , 



(14) 



and the term Eg comes from the expansion of the 
Born cross-section in powers of 0^, 



.(15) 



The remaining contributions to E in (12) are con- 
sidered below. 

2.1. Single hard-photon emission 

In order to calculate the contribution to E due 
to the hard-photon emission we start from the 
corresponding differential cross-section written in 
terms of energy fractions xi,2 and transverse com- 
ponents qi2 of the final particle momenta [^: 



■*e e 7 



dxid^qj'-dx2d^q^ tt 

R{xi]qj;,q^2) ^(1 - X2) 
{qiY {l-U{-(qif )f 

R{x2\q2.qi) ^(1 - xi) 



{qiY {l-Ii{-{qi f ))\ 

2 ^{qif{l-xf 



(16) 



where 

R{x\qi,qi 



1 
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2m2(l - xfx {di - d2) 
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d\dl 



(17) 
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1 + .x^ 
di = m^{\-xf + {qi 
d2 = m^{l-xf + {qi 

and we use the full photon propagator for the t- 
channel photon. Performing a simple azimuthal 
angle integration of Eq. (|l6|) we obtain for the 
hard-photon emission the contribution E^: 



l-A 



a 

TT 



Ax 



l + x' 
l-x 



Az 



z2(i-n(-zg2))2 



X + Q{x^p^ - z)] (L - 1) + k{x, z) 



(1 



k{x, z) = \ ^ [1 + Q{x^p' - z)] + Li 

1 + 

+Q{x^p^ -z)L2 + Q{z - x^p^)L3 , 



(18) 



where L = ln(zQl/ni^) and 
- l)(p2 _ z) 



In 



In 



In 



(a; — z)(xp'^ — z) 
{z — x^){x^p'^ — z) 



x'^{x — z){xp'^ — z) 
{z — x2)(xp2 - z) 



[x — z){x'^p'^ ~ z) ^ 

It is seen from Eq. (|8|) that E^ contains the aux- 
iliary parameter A. This parameter disappears, 
as it should, in the sum YP = -I- S^+^, where 
j^v+s jg j^Yic contribution of virtual and soft real 
photons: 



^ /d2;(i-n(-zg?))-2 

1 

X l^{L-l)P{x)[\ + Q{x^p^ - z)] 
l + x^ 



l-x 



■k{x, z) — (5(1 — a;) 



(20) 



where is the non-singlet splitting kernel. 

3. Radiative corrections to 0{a^) 

A systematic treatment of all 0(0^) contribu- 
tions is absent up to now. This is mainly due to 



the extreme complexity of the analysis (more then 
100 Feynman diagrams are to be taken into ac- 
count considering elastic and inelastic processes). 
Nevertheless in the case of small scattering angles 
we may restrict ourselves by considering only di- 
agrams of the scattering type. 

3.1. Virtual and soft corrections to the 
hard-photon emission 

By evaluating the corrections arising from the 
emission of virtual and real soft photons which 
accompany a single hard-photon we will consider 
two cases. The first case corresponds to the emis- 
sion of the photons by the same fermion. The 
second one occurs when the hard-photon is emit- 
ted by another fermion: 



= da«(^+^)+d(Tff 



H{S+V) 



(S+V) 



, , , (S+V) 



(21) 



In the case when both fcrmions emit, one finds 
that: 



yH ,y{S+V) _ r,yH('^\ 



(L- l)lnA 



(22) 



where T,^ is given in Eq. (|T^. A more complex 
expression arises when the radiative corrections 
are applied to the same fermion line. Here the 
cross-section may be expressed in terms of the 
Compton tensor with an off-shell photon jl^, 
which describes the process 



7*(g) + e (pi) -> e (qi) +jik) + (7,„ft)- 
The result has the form: 



(23) 



l-A 



L < 21nA-ln; 



da;(l + a;2) 
l-x 



X [{L - 1)(1 + 6) + fc(.T, z)] + -\ii^x 
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+ (1 + e)[-2 + Inx - 21n A] + (1 - 9) 



—L In X + 2 In A In a; — In x ln(l 



~Li2(l - x] 

jl-x)^ 
"2(1 + ^2) 



x{l — x) + 4a; In a; 
2(l + a;2) 



In^ X 



(24) 



Li2(a;) = - / jln{l-t), 



where k{x,z) is given in Eq. (|l 
e(a;2p2_2). 



and 8 = 



3.2. Double hard-photon bremsstrahlung 

We now consider the contribution given by the 
process of emission of two hard photons. We will 
distinguish two cases: a) the double simultane- 
ous bremsstrahlung in opposite directions along 
electron and positron momenta, and b) the dou- 
ble bremsstrahlung in the same direction along 
electron or positron momentum. The differential 
cross-section in the first case can be obtained by 
using the factorization property of cross-sections 
within the impact parameter representation [p^ . 
It takes the following form [n3| : 



^^e+e ^(e+7)(e 7) ^ q, 

dxid'^qidx2d'^qi ~ J 7r(fe^)4 
^^ Rixi;qi,k^)R{x2;qi,-k^) 

(i-n(-(fe^)2))2 



(25) 



where R{x\q^.,k ) is given by Eq. (0). The 
calculation of the corresponding contribution 
to S is analogous to the case of the single hard- 
photon emission and the result has the form: 



1 



dz 



z2(i-n(-zg2))2 



l-A 

da;i 



i-A 



dx-. 



1 



1 



l-xi 1 



$(a;i,z)$(x2,z), (26) 



where (see Eq. ([Tg|)): 

$(a;, z) = (L- l)[e(z - l)e(p^ - z) 
+e(z-a;2)e(p\2„2)] 

-^Lt\-Q{x^ - 2) + ©(^ - x^P^)\ 



Li 



1 + X? 
(1 - xf 
1 + a;2 



-H(e(l-z)-e(z-p2))ij^ 



e(z - a;'')e(a;>^ - z) 

e(z-i)e(p2-z) 

(z-a;)(p2 



(xp2 — z){z — 1) 

Let us now turn to the double hard-photon 
emission in the same direction and the hard e"*" 
e~ pair production. Here we use the method de- 
veloped by one of us [|l6|,0. We will distinguish 
the collinear and semi-coUinear kinematics of fi- 
nal particles. In the first case all produced par- 
ticles move in the cones within the polar angles 
Bi < ^ \ centered along the charged-particle 
momenta (final or initial). In the semi-coUinear 
region only one produced particle moves inside 
those cones, while the other moves outside them. 
For the totally inclusive cross-section, such a dis- 
tinction no longer has physical meaning and the 
dependence on the auxiliary parameter Bq disap- 
pears. We underline that in this way all double 
and single-logarithmical contributions may be ex- 
tracted rigorously. The contribution of the region 
when both the photons move outside the small 
cones does not contain any large logarithm L. 
The systematic omission of those contributions in 
the double bremsstrahlung and pair production 
processes is the source of uncertainties of order 
(a/7r)2 < 0.6 • 10-5. 

The contribution of both collinear and semi- 
collinear regions (we consider for definiteness the 
emission of both hard photons along the electron, 
since the contribution of the emission along the 
positron is the same) has the form: 



dz 



1-2A 1-x-A 



dx 



dxi 



z2(l-n(-zQ2))2 



A 



a;i(l — a; — a;i)(l — xiY 



B 



(27) 



I"" = A e(a;2p2 _ ^ 

+CQ{{l-xifp'~z), 

where B and C are known functions [Q. 

The total expression E^t^ which describes the 
contribution to ( [T^ from the two-photon (real 
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and virtual) emission processes reads as follows: 

(28) 



C = IniA. 

The leading contributions YP^ , have the fol- 
lowing forms: 

2^n' J 



X I dx\ -P^^\x) [e{x^p^ - z) + l] 



oc 



(29) 



dxi 



dx2P{xi)P{x2) 



Xc/Xi 



X [e(z - i)e(p2 - z) + e(z - xj)e{xjp^ - z)] 
X [e{z - i)e(p2 - z) + e(z - xDeixlp"^ - z)] . 

We see that the leading contributions to E^''' 
may be expressed in terms of kernels for the evo- 
lution equation for structure functions. 

The functions 0''"'' and (jj^f in expression 
Eq. (^) collect the next-to-leading contributions 
which cannot be obtained by the structure func- 
tions method They have a form that can be 
obtained by comparing the results in the leading 
logarithmic approximation with the logarithmic 
ones given above. 

4. Pair production 

Pair production process in high-energy e"*" e~ 
collisions was considered about 60 years ago (see 
]l3| and references therein). In particular it was 
found that the total cross-section contains cubic 



terms in large logarithm L. These terms come 
from the kinematics when the scattered electron 
and positron move in narrow (with opening an- 
gles ~ m/e) cones and the created pair have the 
invariant mass of the order of m and moves prefer- 
ably along either the electron beam direction or 
the positron one. According to the conditions 
of the LEP detectors, such a kinematics can be 
excluded. In the relevant kinematical region a 
parton-like description could be used giving 
and iy-enhanced terms. 

We accept the LEP 1 conventions whereby an 
event of the Bhabha process is defined as one in 
which the angles of the simultaneously registered 
particles hitting opposite detectors (see Eq. (|4^)). 

The method, developed by one of us (N.P.M.) 
p6| , p^ , of calculating the real hard pair pro- 
duction cross-section within logarithmic accuracy 
(see the discussion in sect. 6) consists in separat- 
ing the contributions of the coUinear and semi- 
coUinear kinematical regions. In the first one 
(CK) we suggest that both electron and positron 
from the created pair go in the narrow cone 
around the direction of one charged particle [the 
projectile (scattered) electron (q^) or the pro- 
jectile (scattered) positron (92)]- 

(r6'o/m > ^ Pi, P2, 91, 92 ■ (30) 

The contribution of the CK contains terms of 
order {aL/n)^, {a/n)^ Lln{0o/e) and ia/Tr)^L, 
where 9 = pTq-^ is the scattering angle. In the 
semi-coUinear region only one of conditions (|30| ) 
on the angles is fulfilled: 



p+p < Oo, P±P, > Oo ; 
or pZpi<Oo, p^i>eo: 

or pZPi > Oo, pTpi < 00 . 



(31) 



The contribution of the SCK contains terms of 
the form: 



iln 



70 



L. 



(32) 



The auxiliary parameter 9o drops out in the total 
sum of the CK and SCK contributions. 

Taking into account the leading and next-to- 
leading terms we can write the full hard pair con- 
tribution including also the pair emission along 
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the positron direction, after the integration over 
X2 as 

l-A 

... = 2-;|/J / d.{L^(l + e)i?(x) 



^1.2 (a;) 
d{x) 

R{x) 



£[eFi(x)+F2(x)] 

d{x)+Ci,2{x), 
1 



(33) 



1 - 
1 1 



3ln(l 



20 



&x 



'-{A + lx + Ax"^) 



3 1 - X 
+ (1 + x) Inx, 

where C\ {x) and C2 (x) are known functions 0] . 

Eq. (|3^ ) describes the small-angle high-energy 
cross-section for the pair production process, pro- 
vided that the created hard pair can move along 
both electron and positron beam directions. 

The contribution to the cross-section of the 
small-angle Bhabha scattering connected with 
the real soft (with energy lower than Ae) and 
virtual pair production can be defined pl| ] by the 
formula: 



C 
20 



1 
17 

InA^ 



3 



■In^ A 



(34) 



Using Eqs. (^ and (|34|) it is easy to verify that 
the auxiliary parameter A is cancelled in the sum 



Chard + CT.oft^ 



We can, therefore, write 



the total contribution cTp^i, as 



dz 



L2(i + _ln(l-a;c) 



-^ln(l - Xc 



^ln^(l -Xc) 



^e.(^-^ln(l-x)) 
1-x ^9 3 ^ " 



dx\L^(i + e)i?(x) + £(eci(x) 



+ C2{X))\ , 



(36) 



R{x) R{x) 



3(1 -x)' 



e = i-e. 



The right-hand side of Eq. (^) gives the con- 
tribution to the small-angle Bhabha scattering 
cross-section for pair production. It is finite and 
can be used for numerical estimations. The lead- 
ing term can be described by the electron struc- 
ture function Dl{x) Jl^ . 

5. Terms of 0{aCf 

In order to evaluate the leading logarithmic 
contribution represented by terms of the type 
(a£)^, we use the iteration up to 0^ of the mas- 
ter equation Q obtained in Ref. ||l^. To sim- 
plify the analytical expressions we adopt here a 
realistic assumption about the smallncss of the 
threshold for the detection of the hard subpro- 
cess energy and neglect terms of the order of: 



x'li-Cf < 3 - 10 



-5 



71 = 1,2,3. 



(37) 



We may, therefore, limit ourselves to consider 
the emission by the initial electron and positron. 
Three photons (virtual and real) contribution to 
E have the form: 



I]3^ = \{^Cfj J y"dxi y"dx2 e(xiX2 - Xc) 



Xr. Xf. 



\b{\~X2)P^^\x^)Q{x\p' -z) 
6 



^^P^''\xi)P{x2)QlQ2 



(l + 0(x3)), (38) 



6162 



e z 



X2 N 



(39) 
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P'^'(a;) is the three-loop spUtting function. The 
contribution to E of the process of pair produc- 
tion accompanied by photon emission when both, 
pair and photons, may be real and virtual has the 
form (with respect to paper by M. Skrzypek 
we include also the non-singlet mechanism of pair 
production): 



E^^^ -y^li^CfJ dzz-^j dxij dX2 

1 Xc Xc 

xQ{xiX2 - Xc) iliR'^ixi) - ^R%Xi)] 



x6il~X2)e{xlp^ - z) 

1 



2 Pix2)Rixi) 6162}, 
ZX-i 



(40) 



where 



Rix) = R%x) + -P{x), 

R'ix) = ^—^{4: + 7x + 4x^) + 2{l + x)lnx, 
3x 

R^{x) = R'{x){^+2lii{l-x)) 
+ (l + x)(-ln2a; + 4Li2(l-a:) 
+ l{-9-3x + 8x^)\nx 

+ ^{---8 + 8x + 3x^) + ^P^^\x). (41) 

6 X 6 

The quantity S depends on the parameters Xc, p 
and Qf. 



6. Estimates of neglected terms and nu- 
merical results 

The uncertainty of our calculations is defined 
by neglected terms. Let us list them. 

a) Terms of the first order RC coming from 
annihilation-type diagrams (15): 



^2 

J ^Ae< 0.10- 10-\ (42) 
el 

b) Similar terms in the second order do not 



exceed (see sect. 4) 

{^?el |^^4<o23.10-^ (43) 



{-Y{e^-e{)c* < 0.5- 10-^ 

TT 

c) We neglect terms which violate the eikonal 
approximation: 



(44) 



^ ^ < 0.3. 10-6. 

TT S 

d) We omit term of the second order which are 
not enhanced by large logarithms: 



{-Y = 0.5 • 10" 



(45) 



e) Creation of heavy pairs {fxfi, tt, tttt, . . .) 
gives in sum at least one order of magnitude 
smaller than the corresponding contribution due 
to light particle production p^ : 

^TTTT + S,,,, + < 0.1 T."^"' < 0.5 • 10"^ (46) 

f) Higher-order corrections, including soft and 
coUinear multi-photon contributions, can be ne- 
glected since they only give contributions of the 
type (aL/nY < 0.2 ■ 10~5 or less. 

g) The terms in the third order associated with 
the emission off the final particles^: 

Xc( — )^ < 0.3- 10-^ (for = 0.5). (47) 

TT 

Regarding all the uncertainties a)-g) and (82) 
as independent ones we conclude the total theo- 
retical uncertainty of our results to be ±0.006%. 

Let us define Eq to be equal to So|n=o (see 
Eq. (21)), which corresponds to the Born cross- 
section obtained by switching off the vacuum po- 
larization contribution H{t). For the experimen- 
tally observable cross-section we obtain: 



Eg (1 + ^0 + ^^+^'^+ -5^^' 



* Usually, in a calorimetric experimental set-up such terms 
do not contribute. 



(48) 
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where 

I][; = So|n=o = l-p"' + Sw' + Se|n=o (49) 
and 

So-^^,S^-§^,S^-'-^,---. (50) 

Zjg Zjq Zjq 

The numerical resuhs are presented in Table 1 . 
Each of the contributions to a has a sign that 
can change because of the interplay between real 
and virtual corrections. The cross-section cor- 
responding to the Born diagrams for producing a 
real particle is always positive, whereas the sign of 
the radiative corrections depends on the order of 
perturbation theory. For the virtual corrections 
at odd orders it is negative, and at even orders it 
is positive. When the aperture of the counters is 
small the compensation between real and virtual 
corrections is not complete. In the limiting case 
of small aperture (p — *■ 1, Xc 1) the virtual 
contributions dominate. 

The numerical results were obtained by using 
the NLLBHA fortran code ||. 

The approach described above to the small 
angle electron-positron cross-section can be also 
used to evaluate, with a next-to-next-to-leading 
accuracy, radiative corrections to the electron line 
at HERA in the small-x region. This problem is 
under investigation. 

The analytical and the numerical calculations 
for the cross-section in the non symmetrical 
Narrow- Wide configuration are in progress and 
will be presented elsewhere. 
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